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We study non-equilibrium transport through a spin-1 Kondo dot in a local magnetic field. To 
this end we perform a two-loop renormalization group analysis in the weak-coupling regime yielding 
analytic results for (i) the renormalized magnetic field and the g-factor, (ii) the time evolution of 
observables and the relevant decay rates, (iii) the magnetization and anisotropy as well as (iv) the 
current and differential conductance in the stationary state. In particular, we find that compared to 
a spin- 1/2 Kondo dot there exist three additional decay rates resulting in an enhanced broadening 
of the logarithmic features observed in stationary quantities. Additionally, we study the effect of 
anisotropic couplings between reservoir and impurity spin. 
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I. INTRODUCTION 

The Kondo model consists of a localized spin- 1/2 which 
is coupled to the electrons in the host metal via a spin ex- 
change interaction. The original motivation to studjdl 
the properties of this and related models was the obser- 
vation of a resistance minimum in some metals in the 
1930s. More than two decades ago it was realizecP^^^ 
that the Kondo model can also be applied to describe 
transport experiments through quantum dots, where the 
localized spin-1/2 models a minimal two-state quantum 
dot which is coupled via exchange interactions to two 
(or more) leads held at different chemical potentials. 
The tremendous developments in the ability to engineer 
nanoscale devices also led to the experimental observa- 
l-jpjpHU! Qf Kondo physics in transport experiments. This 
in turn triggered many theoretical studiea^^HSSof the non- 
equilibrium transport properties of Kondo quantum dots. 

It is well-known that the Kondo model with larger spin 
on the dot show s-- strong deviations from standard 
Fermi liquid behavior at low temperatures. For example, 
at low temperatures the local spin is not fully screened by 
the electrons and a residual magnetization remains. De- 
spite the recent observation of this underscreened Kondo 
effect in transport experiments in single molecule#2rt35^ 
the non-equilibrium transport properties of the spin-1 
Kondo model have not been studied theoreticallySSEl] j-^ 
the same detail as its spin-1/2 counterpart. 

In this article we fill this gap by performing a two-loop 
renormalization group (RG) analysis for a spin-1 Kondo 
dot with finite bias voltage. To this end we apply the real- 
time renormalization group method in frequency spac^^ 
(RTRG-FS) to derive analytic results for the renormal- 
ized magnetic field and g-factor, decay rates and time 
evolution as well as magnetization, anisotropy and dif- 
ferential conductance in the stationary state. All results 
are obtained in a systematic weak-coupling expansion in 
the renormalized exchange coupling J(A) ^ 1. The lat- 
ter condition is satisfied as long as the maximum of all 



occurring energy scales is much larger than the Kondo 
temperature T^, 

Ac = max{y, ho} > Tk, (1) 

where V and are the applied bias voltage and mag- 
netic field, respectively, and Tx denotes the Kondo tem- 
perature at which the system enters the strong-coupling 
regime. We find that the properties of the spin-1 Kondo 
model are similar to the spin-1/2 system with the main 
differences (i) that the larger local Hilbert space on 
the dot results in additional (quintet) decay rates be- 
side the usual longitudinal and transverse spin relaxation 
rates, (ii) that these rates describe the decay of the spin 
anisotropy T oc S''S^ (which is trivial in the spin-1/2 
model), where S"* are the components of the spin opera- 
tor on the dot, and (iii) that the existence of these addi- 
tional rates yields a broadening (relative to the spin-1/2 
model) of logarithmic features in stationary quantities 
like the susceptibilities and the differential conductance. 

The paper is structured as follows. In the next two sec- 
tions we define the spin-1 Kondo dot and briefly review 
the applied RTRG-FS method. In Sec. |IV| we present the 
results for the renormalized magnetic field and g-factor, 
decay rates and time evolution as well as magnetization, 
anisotropy and current in the stationary state. In partic- 
ular, we highlight the differences to the spin-1/2 model 
studied in Ref. [25l In Sec. |V]we investigate the effect of 
anisotropic exchange couplings. Our calculations mainly 
follow the analysis of the spin-1/2 Kondo modeP^ we 
discuss the required modification for the study of the 
spin-1 system in the appendix. 

II. SPIN-1 KONDO DOT 

In this article we consider a spin-1 Kondo dot, i.e. 
a quantum dot whose internal degree of freedom is 
given by a spin S = 1, which is attached to two elec- 
tronic reservoirs via a spin exchange interaction J. This 




FIG. 1. (Color online) Sketch of the spin-1 Kondo model 
([2|. The spin-1 on the quantum dot is coupled via exchange 
interactions J to the electron spins in the two reservoirs at 
different chemical potential ^l/h = ±1^/2. 



set-up is rele vant for recent transport exp eriments on 
Cgo molecule^SlESl ^^^j cobalt complexe^. The spin- 
1 Kondo model can be regarded as a direct general- 
ization of a spin- 1/2 Kondo dot studied previously in 
Ref. 1^ which itself is related to the single-orbital An- 
derson model by a Schrieffer- Wolff transformation. 

Specifically, the Hamiltonian of the spin-1 Kondo dot 
is given by (see Fig. [T]) 



H = '^{ek- fia) ctka^akcr + flQ S*^ 



(2) 



fc'cr' 



aa' kk' era' 



Here c^^^ and Caka create and annihilate electrons with 
momentum k and spin a ^ in lead a — L, R, and 
(T denotes the Pauli matrices. The chemical potentials 
in the leads are given by /i^/i^ — zLV/2 thus applying a 
finite bias voltage V = — jJ-R to the dot. Furthermore 
we introduce an ultra-violet cutoff D via the density of 
states 



N{uj) 



(3) 



and note that the exchange coupling Jo is dimensionless 
in our convention. The spin on the dot is subject to an 
external magnetic field ho applied along the z-direction. 
If we use its direction as the quantization axis, the com- 
ponents of the spin-1 operator on the dot are explicitly 
given by 



S 



The impurity spin on the dot is coupled via an antifer- 
romagnetic exchange interaction Jq > to the electron 
spins in the two reservoirs. (We note that the situation 
of reservoir-dependent exchange couplings Ja/s follows as 
a straightforward generalization as long as the condition 




Jlr = JllJrr is satisfied, 
ture and ho,V >0. 



We consider zero tempera- 



III. RTRG-FS ANALYSIS 

As it is well known for Kondo models like ([2| perturba- 
tion theory in the exchange coupling Jq leads to logarith- 
mic divergencies. In order to takle this problem we use 
the RTRG-F^22l niethod, which is particularly suited to 
derive analytic results for the spin decay rates, the renor- 
malized magnetic field, the dot magnetization as well as 
the differential conductance in non-equilibrium. Our cal- 
culations follow the analysis of the spin- 1/2 Kondo model 
performed in Ref. 1251 we present here a brief summary 
to set up the notations and discuss the essential modifi- 
cations required for the spin-1 system in the appendix. 

We start with the time evolution of the density matrix 
of the full system p(i), which is governed by the von 
Neumann equation 

p{t)^-nH,p{t)]. (5) 
This equation is formally solved by 

p{t) = e-''^(*-*«)p(io)e'^'*"*''' = e-''^(*-*«V(io), (6) 

where the Liouvillian L is defined via L . = [H, . ] . We 
assume that for t < to the reservoirs and quantum dot 
are decoupled and hence that the initial density matrix 
at t — to is of the form p{to) = Psito)PresP?es- Here p"^g 
is the grand canonical density matrix of reservoir a in- 
cluding the chemical potential /1q and psito) denotes the 
initial density matrix of the dot system. The central ob- 
ject of our analysis is the reduced density matrix of the 
dot ps{i)i which is obtained by tracing out the reservoir 
degrees of freedom. Instead of studying the time evolu- 
tion of the density matrix directly we further perform a 
Laplace transform 



ps{z)^ / dte^^(*-*°)ps(t) 



'to 
Tr 



z- L{z) 



Psito)Pres, 



(7) 



where z is the Laplace variable. In order to calcu- 
late ps{z) we expand in the exchange interaction, which 
yields a diagrammatic representation in Liouville space. 
Resumming the resulting series one obtaind^ 



Ps{z) 



Lf{z) 



PS (to), 



where the effective dot Liouvillian is given by 



Lfiz) = L^->+nz)- 



(8) 



(9) 



Here the first term . — Hq [S^, . ] describes the action 
of the applied magnetic field on the spin, whereas the sec- 
ond term S(z) incorporates all interaction effects of the 
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coupling to the leads. From L'g^{z) we obtain the renor- 
malized magnetic field, the decay rates, and the reduced 
density matrix of the dot (and thus the dot magnetiza- 



tion) which are presented in Sees. IV A IV B and IV C 



respectively. A similar analysis can be performed for the 
current kernel yielding the current and differential con- 
ductance through the system (see Sec IV D ) . 



The RTRG-FS allows to calculate the effective dot Li- 
ouvillian and current kernel as well as other observables 
(see e.g. Ref. [26l in a controlled weak-coupling expansion 
in the renormalized ex change couplings J. For Kondo- 
type models this is don^SHHIlll ]-,y fji-g^^ deriving and solv- 
ing the poor-man's scaling (PMS) equations. In the sec- 
ond step one then expands all quantities of interest in 
the effective coupling J(A) while staying in the weak- 
coupling regime J ^ 1. 

In a more physical picture, the KG procedure reduces 
the effective band width by integrating out the high- 
energy degrees of freedom. To accomplish this the ini- 
tial band width D is replaced by an effective band width 
A which flows from Aq ~ D to the physical scale Ac 
which in turn is determined by the restriction to the 
weak-coupling regime. The RG equations describing the 
flow of the interaction vertices and effective Liouvillian 
are obtained from the generic RG equations derived in 
Rcf. "2P by specifying to the spin-1 case (see appendix 
for details). In particular, the renormalization of the ex- 
change coupling in leading order is governed by the PMS 
equation 



dJ 
dA 



A 



J' 



(10) 



which is identical to the spin-1/2 model. The solution 
reads 

•^^^)^^Iir^' TK^Aoe-'/''°, (11) 

where the Kondo temperature is defined as the scale at 
which the effective coupling constant diverges. The scal- 
ing limit is obtained by taking Aq — > oo and Jq — t- while 
keeping constant. Below we will always perform the 
scahng limit unless stated otherwise [e.g. in (15)]. 

As the PMS coupling diverges at A = Tk, the weak- 
coupling condition J(A) <^ 1 can be translated into a 
condition for the physical cut-off scale 



A. 



max{|z|, y, ho}, 



(12) 



which has to satisfy Ac ^ Tx- As z — >■ for station- 
ary quantities we thus recover the condition ([T]). In this 
regime a perturbative expansion of the Liouvillian and 
the current kernel in Jc = J(Ac) is well-defined. Hereby 
the precise value of A^ is irrelevant^^ '^^ as changes 
Ac — > Ac with Ac/ Ac ^ 1 only lead to higher-order cor- 
rections in Jc which will not be analyzed. The explicit 
calculation of the Liouvillian and current kernel follows 
along the lines of Ref . HHI with all necessary modification 
required for the spin-1 Kondo model discussed in App.[B| 



In the following we present the results for the renormal- 
ized magnetic field and g-factor, decay rates and time 
evolution, magnetization and anisotropy as well as the 
differential conductance through the system, where we in 
particular focus on the differences to the spin-1/2 model. 



IV. RESULTS 



The systematic weak-coupling expansions for the ef- 
fective dot Liouvillian Lg^iz) and the current kernel 
I]^(z) involve terms of the order In, where In de- 
notes terms of the form ln[iAc/ (z -I- A -I- iTi)] with some 
physical scale A = V, hp, V ± ho and F^ ^ Jc A the de- 



cay rates (see Sec. IV B I. From the dot Liouvillian one 



can derive the renormalized magnetic field as well as 
the decay rates by determining the poles of the resol- 
vent l/[z - Lf{z)]. If we use the spectral decompo- 
sition L'g^{z) = Ylii Pi{z), where Xi{z) and Pi{z) 
denote the eigenvalues and correspondin g proj ectors re- 
spectively, they are determined by solving^SBH 



(13) 



where the real parts of Zi yield oscillation frequencies in 
the time evolution of observables and thus the renormal- 
ized magnetic field, while the imaginary parts lead to ex- 
ponential decay and thus determine the decay rates (see 
App.[A|for details). Furthermore, the stationary reduced 
density matrix p|* is obtained by solving 



(14) 



thus yielding the stationary dot magnetization and 
anisotropy (see Sec. |IV C I . Finally the current through 
the system is derived from pg and the current kernel. 



A. Renormalized magnetic field and g-factor 

The renormalized magnetic field h, which emerges from 
the externally applied magnetic field ho due to screening 
processes from the electrons in the leads, is determined 
by the real parts of the poles Zi derived from the self- 
consistency equation s (|13[ ). We present explicit relations 
for the ZiS in Eqs. (|Al3l-(|Al8|). To first order in the 



exchange coupling we find 



/iW = [l-(Jc- Jo)]/lo 



(15) 



which is identicaP^I^IIlSl iq |;]^g renormalization of the 
magnetic field in the spin-1/2 Kondo model. Differences 
in the models manifest themselves in the logarithmic cor- 
rections to (15), which yield 



/id 



J' 



J: 



hC-{h) + ^{V ~h)C-{V -h). (16) 



In the derivation of (161 we have neglected all terms in 
order that do not contain logarithms at either h = 



3 




system. As the local Hilbert space as well as the corre- 
sponding Liouville space of the spin-1 dot has a larger 
dimension than its spin- 1/2 analog, we expect the ap- 
pearance of additional decay rates in the spin-1 case. In 
the following we first analyze these rates in the case of 
vanishing magnetic field and then consider the case of a 
finite field. 



1. No magnetic field ho — 



FIG. 2. (Color online) The g-factor g — 2dh/dho for the spin- 
1 (solid line) and the spin-1/2 (dashed line, see Ref.I25| rnodel 



for V = 100 Tff, where Jc — J(Ac) is given in Eq. (Ill with 
Ac = max{\/, h}. We observe a logarithmic feature a,t h — V 
which is less pronounced in the spin-1 model. 



or /i = performed the scaling limit Jq 0, and 
further replaced ho ^ h (which only leads to higher- 
order corrections). In addition, we have introduced the 
short-hand notation 



£_(a;) =ln 



A. 



In 



A, 



n)2 



(17) 

where the appearing decay rates F^, F^, and F?^will be 
discussed in detail in the next section (see Fig. ^ . Here 
we already note that at resonance V = h they satisfy 
Tl ^ rl and Fj 7^ F^, i.e. ^ is well-defined. 



From the renormalized magnetic field ( 16 1 we can di- 
rectly obtain the g-factor 



(V-h), (18) 



which is plotted in Fig. [2] Compared to its bare value 
the g-factor is reduced due to the screening of the spin 
on the dot by the electron spins in the leads. We further 
observe two logarithmic features at ft, = and h = V. 
Due to the reduction of the renormalized magnetic field 
compared to ho the latter appears at ho > V, however, 
up to logarithmic corrections in order the position 
equals the spin-1/2 situation. Differences show up in the 
line shapes close to ft. = and h = V, where the spin- 
1 result clearly shows a larger broadening. The reason 
is that in the spin-1 case both broadenings Fj and Fj 
(which satisfy Fj > F^), appear in the logarithms, while 
in the spin-1/2 case one findPI rl = T]. Away from 
these positions, however, the g-factors are identical. 



B. Decay rates and time evolution 



For vanishing magnetic field Hq = the model is fully 
spin isotropic. Thus the local Liouville space of the dot 
can be decomposed into irreducible representations of 
SU(2). For a spin-1/2 system these are a singlet and 
a triplet (according to 5<Xi^ = 0®1), while for the 
spin-1 system one finds a singlet, triplet, and quintet 
(l(g)l = 0®1©2). Therefore the effective dot Liou- 
villian of the spin-1 model ^ decomposes as 

Lf{z) = -ifsizW ~ ^It{z)L' - if,{z)L\ (19) 

where the superoperators , L*, and L'' project onto 
the respective subspaces (see App. |A]). In leading order 
the coefficients in (19) read fs{z) = 0, ftiz) = ttJ^V, 
and fq{z) = SttJ^V. As there is no dependence on the 



Laplace variable z the self-consistency equations ( 13 1 are 

s, t, q) with 



trivially solved by Zi 



-iVi {i 



F« 



0, F, = ttJ^V, 



SttJ^V. 



(20) 



Obviously the corresponding subspaces in Liouville space 
are 1-dimensional (F^), 3-dimensional (Ft), and 5- 
dimensional (Fg). The difference to the spin-1/2 modeP^ 
is the appearance of the additional rate F^. Similarly the 
decay rates in a general spin-S' Kondo dot in the absence 
of a magnetic field are given by 



l)JlV for J =0,1, 



,2S. 



(21) 



In order to interpret the result ( 20 ) we investigate the 



time evolution of the reduced density matrix, which is 
given by the inverse Laplace transform of ([t]) 



PD{t) = 



2tt 



dz e 



-oo+io+ z — L'j^{z) 



PDiO). (22) 



The initial density matrix is assumed to be of the form 

p{to = 0) = Ps{0)PresP?es ^ith 



Pll PW Pl-1 

Ps{0) = I Pqi Poo Pq-i 
P-ii P-io /O-i-i 



(23) 



In this section we investigate the imaginary parts of 
the poles Zi obtained by solving (13 1. As we will show 



they lead to an exponential decay in the time evolution 
of observables and thus constitute the decay rates of the 



We perform the integral in ( 22 ) by using the decomposi 



tion ( 19 ) for the Liouvillian. In leading order this yields 



Ps{t)^[L'+e-^*'L' + e-^^*L'^]psiO). (24) 



4 



We note that the time evolution is governed by purely 
exponential decays. This is due to the absence of a z- 
dependence of L^g in leading order; in general the branch 

cuts of L^g due to logarithmic terms will lead to addi- 
tional power-law corrections^. The time evolution of the 



spin operator is now readily obtained from (24 1 



cx e-^'*, e.g. 
(5^(t)) = (pn-P-i-i)e 



(25) 



Thus the magnetization decays only with the triplet rate 
Ft, which is identicaP^to the spin-1/2 model. However, 
in the spin-1 model one can construct a second non-trivial 
operator on the dot, namely the anisotropy operatoif^S (or 
spherical quadrupole tensor) with components 



1 



2±2 



2±1 



- S^S'^ ± i(S'^5^ - 
-F5^S'^±i(S'^5'^ 



?20 — 



1. 



S-^S-^)], (26) 
^S'Sy)l (27) 

(28) 



(T is relatecP^ to the traceless quadrupole tensor Qij 
introduced in Refs. |46] and |47] to study the trans- 
port through spin valves.) Physically measures the 
anisotropy along the z-direction. (r2o) > shows a ten- 
dency to align the impurity-spin along the z-axis, while 
(^20) < indicates an alignment in the xy-plane. We 
further note that in the spin-1/2 case one finds T20 oc 1, 
i.e. the anisotropy operator is trivial. From (24) we can 
easily infer the time evolution of T 



(r(t)) 



e.g. 



1 



(J2o(t)) = ^ 2poo + P-i-i)e 



-r„t 



(29) 



Thus the situation in the absence of a magnetic field is 
quite simple: The magnetization decays with the triplet 
rate Fj, whereas the quintet rate Fg governs the time 
evolution of the anisotropy operator. In the next section 
we consider how this changes if the spin rotational in- 
variance is broken by the application of a finite magnetic 
field. 



2. Finite magnetic field /lo > 

The breaking of the spin rotational invariance results 
in a decomposition of the triplet and quintet subspaces in 
the local Liouville space of the dot. Hence the dot Liou- 



1000 




FIG. 3. (Color online) The decay rates (|30| for V = 10'' Ta-, 
where Jc = J(Ac) is given in Eq. (Ill with Ac = max{l/, /i}. 
We observe the splitting of the triplet and quintet rate into 
five rates, of which two (Fj and F^) merge for ft > 1/. Inset: 
For /io 'F there are two rates governing the decay along 
and perpendicular to the z-direction. 



respectively) exist, we find five decay rates in the spin-1 
system. As each rate corresponds to a relaxation or deco- 
herence time, the time evolution towards the stationary 
state will be more complex in the spin-1 model (see be- 
low). 

In the small-field limit ^ V the splitting of the 
rates is explicitly given by 

T\^^JlV ^^Jlh, T\=^JlV+ljlh, (31) 



and 



f;; = i^jiv 



irJ^h. 



2 " ' (32) 



We note that the magnetic-field dependence of the triplet 
rates (31) is identical to the spin-1/2 modeP^ with the 
identifications F° — >■ Fi and F^ — > F2. For arbitrary mag- 
netic fields the rates have been obtained by numerically 



solving the self-consistency equations (C1)-(C3|, the re- 
sult is shown in Fig. [3j We observe in particular that two 
of the rates (F" and F^) merge a,t h = V. We note that 



( [Tt] ) is well-defined at resonance as Fj ^ F^ and F J 7^ F^ . 

The time evolution of the reduced density matrix fol- 
lows from ( 22 1 using the spectral decomposition L^g^{z) = 
\i{z)Pi[z) of the dot Liouvilian. If we further approx- 
imate this by expanding around the poles of the resolvent. 



villian does no longer take the simple form ( 19 ) and the 

system possesses more than two decay rates. Specifically Lf{z) « h{zi)Pi{zi) and use Eq. ([l3]), we obtain 
we find that the rates ( 20 ) split up as 



Ft^FO,Fi, F, 



T^O T^l T^2 

^ q-) ^ q-) ^ q- 



(30) 



Ps(t)-5I' 



P^{Zi)ps{t^O). 



(33) 



Compared to the spin-1/2 model, where at finite mag- 
netic field two different rates (the longitudinal and trans- 
verse spin relaxation rates usually denoted by Fi and F2 



Here the real parts of the Zi yield oscillations with the os- 
cillation frequency determined by the renormalized mag- 
netic field, while the imaginary parts correspond to the 



5 



decay rates Q [see Eqs. ( |A13[ )"( |AT8| ]. The first pole 
zi = yields a stationary contribution, i.e. -Pi(O) is the 
projector onto the stationary density matrix ptf- From 



(331 we obtain for example 



(34) 



Corrections to ( pS]) c an be calculated along the lines of 
Refs. EH ESI and 1491 We expect additional terms oscil- 
lating with the freqencies V, h±V, and 2h±V as well 
as power-law corrections to the exponential decays. 

A better understanding of the time evolution and the 
relevant decay rates can be obtained by studying the 
large- field limit, hg ^ V. In this regime we find only 
two effective rates which are related to the five original 
ones by (see inset in Fig. [3]) 

"2 > = 27r 



^ tJ^ q 



-pO pO -pl _ 
^ ^ qT^ q 

[We note that although Fj 



AnJ^h. 



(35) 
(36) 



F^ the logarithmic terms in 



( 16 1 are well defined as we are off resonance.] The time 

(37) 



evolution of the magnetization is now given by 



[Pii - P-i-i 



72 



(e-r^* + e-r^*)(p^io-l-po±i), (38) 



with = 5*^ ± iS'^ and the initial density matrix of the 
dot given by (23), while the anisotropy operator relaxes 

as 



{T2,{t)) 

{T2±l{t)) 
{T2±2{t)) 



1 

71 



(1 



(pii - 2poo + P-i-i 



(39) 

-r^* + e-r^*)(p^io-po±i), (40) 

(41) 



2^2 ^ 



(36 



We observe that the decay of the expectation values of the 
diagonal operators and T20 is purely governed by the 
rate F^, while {T2±2{t)) decays purely with Fj^. Thus we 
conclude that F^ corresponds to the longitudinal spin re- 
laxation rate of the spin- 1/2 system (usually denoted by 
Fi), while Fj^ corresponds to the transverse one (usually 
denoted by F2). Comparing the explicit values (35) and 
with the spin-1/2 result!^ shows that both rates Fj^ 
F^ are twice as large as their spin-1/2 counterparts 
(i.e. their ratio is identical in the spin-1/2 and spin-1 
model). We stress, however, that the time evolution in 
the spin-1 model is more complex, e.g. in the decay of 
(S'±(t)) and (T2±i(i)) both rates F^ and F_l appear. 



C. Magnetization and Anisotropy 

In this section we analyze the stationary values of mag- 
netization and anisotropy before turning to the differ- 
ential conductance in the next section. The stationary 



density matrix of the dot is obtained by solving ( 14 ). In 



troducing the observables magnetization and anisotropy 
by 

M=(5^),t, A=(r2o),t, (42) 
the stationary density matrix can be written as 



M 



20- 



(43) 



The knowledge of the effective dot Liouvillian enables 
us to derive analytic results for the magnetization and 
anisotropy including the leading logarithmic corrections, 
which are given by 



M 



4/1/2 



{h? + i{f2?" 



A 



1 



4(/i)^ 



(44) 



/6(/i)2+3(/2)2 
with 

/i = 2T:Jlh + 2T:JlhCi{h) - 7rj|(F - h)Ci{V - h) 
h = -\jl{2V + Qh+\V-h\,) 

"'\V -h\^-2{V -h)\Ci{V -h) 



-^JlhC,{h)--Jl 



as well as 

Ci{x) 



In 



A, 



in? 



2 I X 
X arctan — r 



In 



arctan 



A,. 



(45) 



(46) 



In the derivation of (44) we have neglected all terms in 
order that do not contain logarithms at either h — 0, 

V = 0, OT h = V. We note that the decay rates Fj and 
F^ cut-off the logarithmic divergencies present in bare 
perturbation theory. As we will see below these two 
rates also appear in the logarithmic corrections to the 
differential conductance. 

Specifically we find that for large magnetic fields h > V 
the spin is completely aligned along the field, i.e. M = 
— 1 and A = 1/a/6. Close to the resonance, i.e. a.t h <V, 

V — h h wc find by expanding (44) up to O(Jln) 



M 



A^ 



1 V 

1~ 



JcV-h 



V6 



1 



h 4 h 
3V-h 3JcV -h 
4 h~ 



Ci{V~h), 



Ci{V~h) 



(47) 
.(48) 



Similarly in the limit of small magnetic fields, ft. <C V^, we 
obtain 



M = 
A = 



8/1 
~3y 
16 



[l + JcCi{h)\, (49) 



[l + 2J,Ci{h)\. (50) 



The magnetization (together with the spin-1/2 result) 
and anisotropy are shown in Fig. [4j 



6 




0.5 1 1.5 2 0.2 0.4 0.6 0.8 1 1.2 

hg/V hg/V 



FIG. 4. (Color online) The relative magnetization Af/Mmin 
for the spin-1 (solid line) and the spin-1/2 (dashed line)'^ 



model for V — 10 Tk, where Jc = J (Ac) is given in Eq. (Ill 
with Ac 
model. 



max{V, /i}. Inset: Anisotropy A for the spin-1 



The corresponding magnetic susceptibility is readily 
obtained from M. At resonance and for small fields we 
find 



DM 
dho 
dM 
oho 



1 



'l + JcCiiV-h)],{51) 



3V 



[1 + Jc£i(/i)]. 



(52) 



FIG. 5. (Color online) Left axis: The susceptibility x = 
for the spin-1 (solid lines) and spin-1/2 (dashed lines)p^ mod- 
els for V = 100 Tk, where Jc = J(Ac) is given in Eq. ( (Tl| ) 
with Ac = ma.x{V,h}. Right axis: where x'"^ is the 

result in 0(Jc), i.e. without logarithmic corrections. We note 
that the logarithmic corrections at ft = and h = V are less 
pronounced in the spin-1 model. 



the following, 1 = 1^. Within the RTRG-FS formalism 
the stationary current is given bj!^ 



(53) 



Except for the logarithmic broadening the susceptibility 
at = y is identical to the spin-1/2 cas^^, whereas in 
the limit h ^ they differ by a relative factor 8/3. We 
note that the same factor also appears in the susceptibil- 
ity of an isolated spin S at finite temperature. In Fig. [5] 
we plot the susceptibility for the spin-1 model as well as 
the spin-1/2 systerrPSl. In order to enhance the visibil- 
ity of the logarithmic terms we further plot x/x^^\ where 
is the result in 0{J^). In comparison to the spin-1/2 
model the logarithmic corrections at — > and h V = 2ttJ^V 



where the current kernel 'S^ (z) can be derived in complete 
analogy to the effective dot Liouvillian. We calculate 
up to order In, thus together with the knowlegde 
of pg up to this order we are able to derive analytic 
results for current and conductance including the leading 
logarithmic corrections. The current is explicitly given by 



I = f[ + MfL + ^Afi. 



where M and A are given by (44) and 



are broadened as they contain two decay rates [see (45 1] 
of which one is larger than the rate appearing in the spin- 
1 /2 susceptibility. We expect this to be a generic feature 
of the spin-S* Kondo dot, i.e. the logarithmic corrections 
will be most pronounced in the 5* = 1/2 model and be- 
come more and more broadened in the large- S' limit. 



D. Differential Conductance 

Finally, we discuss the current and conductance 
through the system. We define the current as the 
change of the electron number in the left reservoir, — 
— ^A^res- Because the number of electrons on the dot 
is fixed, the current in the right reservoir is accordingly 
given by = —I^. Thus we will omit the superscript in 



(54) 



(55) 



fi 



2 I 



{V -h)Ci{V -h) + VCo{V) 
2h~ \V ~h\j) 



(56) 



WCoiV) + AhCiQi) -\V - h\^Ci{V - h) 



fi - ^J! {V - h)C,{V -h)- 2VCa[V) 



with the logarithmic terms given by Eq. ( 45 ) and 
A, 



Co{x) = In 



In 



A, 



(r?)2 



(57) 



(58) 



Fig. |6] shows the differential conductance as a function of 
the applied bias voltage V derived from Eq. ( 54 ) . 
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FIG. 6. (Color online) The differential conductance G = 
in the spin-1 (solid line) and spin- 1/2 (dashed hnej^ model 
for ho — 100 Tk, where Jc = J (Ac) is given in Eq. ( |11[ ) with 
Ac = max{V^,/i}. We observe that the linear conductance 
Giy — >■ 0) in the spin-1 model takes four times the value of 
the one in the spin- 1/2 model [see Eq. ( 59 1] and that the onset 
of inelastic cotunneling processes at V = /i is sharper in the 
spin- 1/2 model. 



We find for small bias voltage F <C /i that there are no 
logarithmic corrections and the current is governed solely 
by elastic cotunneling processes, 



dl 
dV 



(59) 



This is four times the result in the spin- 1/2 model due 
to the additionally available transport channels. (In the 
spin-S dot the elastic cotunneling current will be given 
by /= l{2S+l)S{S + l)J'^V.) 

On the other hand, close to the resonance V k, h we 
observe a jump in the differential conductance which is 
due to the onset of inelastic cotunneling processes (see 
also Refs. [33 EH ES and HO]). Explicitly we find 



G_ 

Go 



27r2 J2 -I- 2Tr^J^Ci{V - h) for V < h, 



^t:^J^ +'^n'J'^Ci{V -h) for V > h, 



(60) 



where we assumed \V — h\ ^ h in both cases, and Gq = 
l/27r is the conductance quantum. Furthermore, from 
(601 we deduce the height of the jump at V — h 



AG 

Go 



'-yj^[i + jcCi{v~h)]. 



(61) 



We note that the analog prefactor in the spin-1/2 model 



equals^ 3/2. The physical origin of this jump is the on- 
set of inelastic cotunneling processes: For finite magnetic 
field the spin states on the dot are separated by the en- 
ergy h. Thus for V < h only elastic cotunneling processes 
are possible in which electrons are transferred from the 
left to the right lead via virtual states on the dot without 



flipping its spin state. Hence the initial and final state of 
the dot have the same energy. For V > h, however, also 
inelastic cotunneling processes will contribute in which 
the spin on the dot is flipped (thus requiring the energy 
h). In analogy to the susceptibility we observe that the 
logarithmic corrections are less pronounced in the spin-1 
model. 

Finally we would like to add that the current noise in 
the spin-1/2 and spin-1 models are qualitatively similar, 
i.e. differences appear only in the decay rates broadening 
the logarithmic corrections^. 



V. ANISOTROPIC MODEL 

So far we have restricted ourselves to the case of 
isotropic exchange couplings. Dropping this restriction 
we have to study the spin-1 Kondo model ^ with the 
generalized exchange interaction 



2 E [4is^^:,'+sycTi,,)+Jss^T^.M 

cxa' kk' aa' 

(62) 

As explained in detail in Refs. |2S] and [H the RTRG-FS 
method integrates out the reservoir degrees of freedom 
in a two-step procedure. In the first step one removes 
the symmetric part of the Fermi function in the reservoir 
contractions in a single (discrete) RG step, which yields a 
well-defined perturbative expansion in Jq and 1/-D, where 
D denotes the band width. For the anisotropic spin-1/2 
Kondo model this step only yields negligible perturbative 
corrections to the Liouvillian and current kernel. In the 
second (continuous) step one then resums the logarith- 
mic divergencies by integrating out infinitesimal energy 
shells in the remaining asymmetric part of the reservoir 
contractions. This yields the PMS equations ( 10 1 as well 



as similar RG equations for the Liouvillian and the cur- 
rent kernel. 

However, the situation is completely different for the 
anisotropic spin-1 Kondo model (62). The presence of a 
finite anisotropy 



{J' 



produces the term [see Eq. (74) in Ref. dS] 



(63) 



(64) 



in the Liouvillian. Here La denotes the superoperator 
corresponding to the square of the local spin. La . = 
[S'^S'^, .]. We note that for spin-1/2 this operator van- 
ishes since S^S^ — 1/4, thus the term (64) does not 
appear. 



Since Eq. (64) contains the initial band width D ^ 
ho,V it will completely dominate the physics of the sys- 
tem. We distinguish two cases: (i) For > the term 



(64) is negative, thus locking the dot in one of the fully 
polarized states with M = ±1 and A = l/-\/6. As the 
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energy spacing to the other states on the dot is propor- 
tional to D only elastic cotunneling processes are possible 
and the current is given by / = tt{Jq)^V. (ii) For < 
the term ( 64 ) is positive and the dot magnetization takes 
the value M = 0. Perturbation theory in the exchange 
couplings then yields / = 0. 

This behavior was previously observed in the 
anisotr opic s pin-1 Kondo model at equilibrium via a PMS 
analysi^il^ as well as an NRG calculatioiP^. In fact, the 
term S^S^ was shown to be a relevant pertur bation on 
the grounds of general symmetry considerationd^i^. 



VI. CONCLUSIONS 

In conclusion, we have studied non-equilibrium trans- 
port through a quantum dot modeled by a spin-1 Kondo 
model in a magnetic field. We obtained analytic results 
for the renormalized magnetic field and g-factor, the mag- 
netization and anisotropy as well as the differential con- 
ductance. The latter shows non-monotonic behavior as a 
function of the bias voltage with a pronounced jump at 
V = h due to the onset of inelastic cotunneling processes. 
In particular, this jump is strongly enhanced by logarith- 
mic corrections as was previously observed in transport 
experiments on carbon nanotube quantum dot^^. 

The transport properties of the spin- 1/2 Kondo dot 
have been studied previously using the same formalisrrP^. 
The main difference in the spin-1 model is the appear- 
ance of three additional decay rates which govern the 
relaxation of the spin anisotropy oc S^S^ on the dot, a 
quantity which becomes constant in the spin-1/2 model. 
Hence the relaxation dynamics of the magnetization and 
current turns out to be much richer in the spin-1 model. 
Furthermore the decay rates appear as cut-off parameters 
in the logarithmic terms of the perturbative expansions. 
As a consequence, the presence of additional rates leads 
to an enhanced broadening of the logarithmic corrections 
of stationary quantities like the susceptibility (see Fig.[5| 
and the differential conductance (see Fig. [6]), which are 
thus less pronounced than in the spin-1/2 model. 

Finally we note that as our approach is based on an 
expansion in the renormalized exchange coupling it is 
not possible to derive reliable results in the low-energy 
regime A^ = maxjl^, h} < Tx where the differences be- 
tween the fully screened spin-1/2 model and the under- 
screened spin-1 model are most pronounced. As the un- 
derscreened Kondo model can be mapped to an effective 
ferromagnetic Kondo model, a perturbative study of the 
transport properties in this regime might be possible us- 
ing a Majorana diagramatic theory as recentlyperformed 
for the investigation of the magnetotransport^. 
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Appendix A: Liouvillian 

Here we present the main modifications required to 
generalize the calculations of Ref. .25, to the spin-1 Kondo 
model. Further details can be found in Ref. [56l 

Superoperators. — The basis for the local Liouville space 
of the dot can be built up by the superoperators and 
defined by 



L+A = SA, L-A = -AS. 
An explicit representation is given by 

L+ = S(g)l, L~=-Cl(g>S 



(Al) 



(A2) 



where C = (1,-1,1). In general the action of a superop- 
erator O on an arbitrary operator B (acting itself in the 
Hilbert space) can be written in a matrix notation as 



{O gs' Ogg' Bgg/ 



(A3) 



Symmetries. — The effective dot Liouvillian obeys the 
same symmetries as the Hamiltonian: hermiticity rela- 
tion, spin-conservation, and probability conservation 

H = ^ [Lf{z)y ^~Lf{-z*), (A4) 

s-s^s-s' ^ (if) = 0, (A5) 

V / ss' ,ss' 

Tr5 (ps(t) - Ps{to)) =O^Tts Lf{z) = 0. (A6) 
Here the c-transform is defined by 

(^'^)ss',sS' — ^*s's,s's i^"^) 

and the trace over the system states is given by 

{TvsL)_^_^,=Y.^--rs-s'- (A8) 

If we choose an ordered basis in the Liouville space of the 
dot by 

{|11) , |00) , 1-1-1) , |10) , |0-1) , |01) , 1-10) , |1-1) , 1-11)}, 

^ (A9) 

the relations ( A5 ) leads to a decomposition of L'g^ into 
one 3x3 block with |s — s'| = 0, two 2x2 blocks with 
\s — s'l = 1, and two 1x1 blocks with \s — s'\ = 2. The 
restrictions ( A4 ) and ( A6 ) then yield further relation^^ 



between the matrix elements in these blocks. 

For vanishing magnetic field /iq = the model is fully 
rotational invariant and the dot Liouvillian decomposes 
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as (19). The projectors on the different spin-subspaces 



are exphcitly given by 



^l+\{L+-L-)\ 



L* = l 



3 ■ 3 
3^ + ^^^"-^^ 
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(AlO) 
(All) 
(A12) 



Poles of the resolvent. — The poles of the resolvent 

S 



l/[z — L'g^{z)] can be parametrized by 



Zl 


0, 




(A13) 


Z2 = 






(AM) 


Z3 = 


-4 = 




(A15) 


2^4/6 = 


— 7* = 




(A16) 


^5/7 = 


~ ^7/5 = 


±h - iV], 


(A17) 


Z8/9 — 


~ ^9/8 = 


±2h ~ iVl, 


(A18) 



where the real parts are related to the renormalized mag- 
netic field h discussed in Sec. |IV A| while the imaginary 



parts yield the decay rates discussed in Sec. IV B 



Appendix B: Evaluation of RG equations 

The derivation of the Liouvillian and current ker- 
nel for the spin-1 Kondo model starts with the results 
(184)-(187) of Ref. [53 For the explicit evaluation we 
decompose the effective dot Liouvillian as L^^(z) = 
Ai(2;)Pi(z), where Xi{z) denote the eigenvalues and 
the corresponding projectors Pi(z)'s are determined by 



the right and left eigenvector^^. The decomposition 
yields for the resolvents appearing in the perturbative 
expansions 



Lf{z) 



1 



\,iz 



■nz) 



1 



-P^{Z^), 



(Bl) 

where we have expanded the eigenvalues and eigenvectors 
around the poles z ^ z,. Eqs. (184)-(187) of Ref. [55] 
can now be evaluated by inserting (Bl) and performing 



the necessary algebra in Liouville space. For example in 
O(J^) logarithmic terms of the form 



^E' 

1=1 



Zi) In 



A. 



'i{z - Zi) 



P^{Z,)G 



(B2) 



appear, where G cc Jc schematically denotes the super- 
operator corresponding to the exchange interaction. We 
proceed by evaluating each te rm in ( |B2[ ) separately. The 
first term with i — 1 vanishe^^siili ^j^g Pi(zi)G — 0. 
This implies in particular that all remaining logarith- 
mic terms are cut-off by a finite decay rate [see ( |A14 )- 
( |A18 )] and thus the perturbative expansion of L'g^ is well- 
defined as long as Jc ^ 1. The terms with z = 8, 9 can be 
calculated straightforwardly. The remaining terms are 
evaluated using the following approximation (which is 
not necessary in the spin- 1/2 model): Introducing the 
short-hand notation 



Kiiz) ^ {z ~ Zi) In 



Zi 



(B3) 



we have to calculate for j = 2, 3 (the remaining pairs 
i — 4:,5 and i = 6,7 can be treated in the same way) 



G 



[H2{z)P2{z2) + H3(z)P3(^3)}G = \G[n2{z) [P2{z2) + ^3(^2)] + U^iz) [^2(^3) + ^3(^3)] ]G 

+ ^g{H2(z) [P2{z2) - P3(22)] - n^iz) [P2{Z3) ' ^3(^3)] } 



G 



\ [n2{z) + n^iz)] G r'^^^^ g + i [^(z) - ^3(^)1 g [P2{z2) - ^3(^2)] g + o{Jl), (B4) 



where we have used G[P2{z) + P3(z)]G = 
GY,l=iPr{z)G = Gl^^^G with l^"^^ denoting the 
projector onto the 3x3 subspace with |s — s'| =0, and 
P^{z■i) = P,{z2) + 0{Jc). Now the first term in ([bI]) 



yields the logarithmic corrections in 0{J^), e.g. the real 
parts of the analog terms with i = 4, 5 and i ^ 6,7 yield 
the logarithmic corrections to the renormalized magnetic 



field ( 16 ) 



On the other hand, the second term in (B4) does not 
lead to logarithmic corrections in 0{J^) as 



H2(z)-H3(z)«ln|| -0(1). 
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(B5) 



The approximation ( |B4[ ) is the main technical modifi- 
cation required for the treatment of the spin-1 Kondo 
model. The possibility to neglect the contributions from 
( B5 ) is essential to derive the analytic results presented in 



the main text. The approximation (B4) obviously yields 
the combination "^2(2) +7^3(2) (simil 
i = 4, 5 and i = 6,7). 



ar for the other pairs 
We note, however, the direct re- 
lations between the functions Hiiz) (and thus the rates 
Ti) and the projectors Pi{zi) are lost. 



We st ress again that our calculations rely on (B5|. 
Thus if ( [B4 ) were to contain 'Hi{z) (as may be the case 
in the calculation of two-point correlation functions) the 
terms of the form ( B5 ) become divergent and the approx- 
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imation ( B4 1 is no longer applicable. 



Finally we note that the imaginary part of ( B4 ) [orig- 



inating from the imaginary part of 'Hi{z)] does not con- 
tain any logarithmic corrections. Thus the logarithmic 
corrections to the imaginary part of L'g^{z), which yield 
the logarithmic corrections to the decay rates and sta- 
tionary quantities, only appear in O(G^). They can be 
evaluated by using a relation similar to (|B4|. 



Appendix C: Self-consistency equations for decay 
rates 



From the poles Zi of the resolvent [see (A14)-(A18|] 



we obtain the following self-consistency equations for the 
decay rates (see Fig. [s]) 



r?/, = 2Ci T V{Ci+C2){Ci-C2), (CI) 
Tl/^ = Cs + lCi + C,Tl v/(C5 + 2C3)2-3(C2)2, (C2) 

(C3) 



= 4^3 + ^4, 



where the right-hand side (RHS) of (CI) has to be eval- 
uated at z = Z2/3, the RHS of (C2) at 2; = Z4/5, and the 
RHS of (C3) at z = zg, respectively, and the appearing 



expressions read 

Ci(z) = '^J^(2\z-h\^+2\z + h\^ + \z-h + V\^ 

+ \z + h + V\^ + \z- h-V\^ + \z + h-V\^y 

C2 ^ 2TTj^h, 

Csiz) = y (21^ - 2/i|2 + \z + V - 2h\^ + \z-V- 2h\, 
Ciiz) 



Jt\2\z ~h\^ + \z^Y -h\^^\z-V - 



^j2(2|z|„ + |z + y|o + |z-F|o). 

The absolute values are broadened by the decay rates 
and explicitly given by ( 46 1 as well as 



2 ( X 
—X arctan —^r 

^ \ r? 
2 



arctan 



-X arctan — 



(C4) 
(C5) 
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